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PREFACE 
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Laboratory  (AFAMRL) ,  Human  Engineering  Division,  Optical  Countermeasure 
program.  This  work  was  performed  under  Contract  F33615-79-C-0500.  The 
Contract  Monitor  was  Mr.  Donald  McKechnie  and  the  Program  Manager  was 
Maj.  Allan  M.  Dickson.  The  SRL  Project  Manager  was  Mr.  Kaile  Bishop. 

The  author  extends  his  appreciation  to  Dr.  Carroll  Day,  Messrs.  Walter 
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Force  Aerospace  Medical  Research  Laboratory,  Wright-Patterson  AFB,  for 
their  valuable  comments  and  discussions. 
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Section  I 
INTRODUCTION 


The  modeling  of  human  performance  in  an  antiaircraft  artillery  (AAA)  system 
has  been  extensively  studied  by  many  investigators  in  the  past  decade  [e.g., 
Kleinman  and  Perkins  (1974),  Phatek  et  al.  (1976),  Kou  et  al.  (1978)],  Most 
of  these  works  dealt  with  the  modeling  of  human  response  in  a  simple  tracking 
task.  In  the  event  of  interrupted  observations,  the  operator's  tracking 
performance  degrades  significantly  during  the  interruption  period  and  poses 
an  appealing  modeling  problem.  The  author  tackled  this  problem  by  degrading 
•several  observer  and  controller  gains  in  the  model  and  proved  to  be  rather 
successful  (Yu  et  al.,  1980).  Efforts  were  then  directed  to  study  human 
response  in  a  manual  tracking  and  firing  task.  In  this  task,  the  operator 
(gunner)  directly  controls  the  gun  turret  and  fires  tracer  rounds  continu¬ 
ously  toward  the  target.  The  gunner  perceives  the  tracer  ending  position 
and  continuously  adjusts  weapon  pointing  in  azimuth  and  elevation  to 
minimize  the  tracer-to-target  error.  In  this  system  mode,  the  gunner  has  to 
play  both  the  role  of  a  tracker  and  a  lead  angle  computer.  The  conventional 
tracking  task  is  greatly  complicated  by  the  inclusion  of  lead  angle  estima¬ 
tion.  Wei  (1981)  developed  an  observer  gunner  model  which  treated  the 
tracer  information  as  delayed  measurements.  The  intent  of  this  paper  is  to 
extend  the  author's  previous  work  to  consider  an  even  more  general  tracking 
and  firing  scenario,  i.e.  to  consider  a  tracking  and  firing  task  subject  to 
external  measurement  interruptions. 

The  interruptions  occur,  in  the  real  world,  through  various  electronic/ 
optical  countermeasures,  weather,  or  terrain  conditions.  In  this  study, 
extensive  manned -simulation  experiments  were  conducted  at  the  Air  Force 
Aerospace  Medical  Research  Laboratory  of  Wright-Patterson  AFB,  Ohio.  A 
typical  helicopter  operational  trajectory  was  used  in  the  experiment.  The 
trajectory  consists  of  three  phases.  During  the  first  phase,  the  target  is 
standing  still  at  certain  altitude  and  is  half  masked  by  some  terrain  con¬ 
figuration.  At  the  onset  of  the  second  phase,  the  target  pops  up  for  a  full 
unmask  flight  and  moves  horizontally.  Blanking  of  target  is  administered  in 
this  phase  only.  Blanking  durations  range  from  1.5  sec,  3.0  sec,  6.0  sec, 
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and  full  blanking.  Certain  repetition  of  blanking  duration  is  also 
included. 

The  structure  of  the  gunner  model  in  a  tracer-directed  fire  system  in  Wei 
(1981)  is  adopted  and  generalized  here.  Nonlinear  ballistic  equation  is 
used  to  compute  the  loci  of  elevation  projectiles.  The  model  consists  of  a 
reduced-order  observer,  a  linear  feedback  controller,  and  a  noise  remnant 
element.  The  remnant  function  which  lumps  all  of  the  random  effects  due  to 
measurement  noise  and  human  neuromotor  response  noise  is  assumed  to  be 
Gaussian  with  its  covariance  being  a  function  of  estimated  target  velocity 
and  acceleration. 

The  effect  of  observation  interruption  is  modeled  by  exponentially  degrading 
the  observer  gain,  the  controller  gains  pertained  to  observed  states,  and 
the  bias  term  in  covariance  function.  Model  parameters  and  time  constants 
are  identified  separately  with  respect  to  no-blanking  and  blanking  empirical 
data  via  a  least-squares  minimization  algorithm.  The  computer  simulation  of 
the  designed  model  shows  that  the  model  predicted  tracking  and  tracer  errors 
are  in  good  agreement  with  empirical  data  over  various  blanking  conditions. 


Section  II 

AAA  TRACER-DIRECTED  FIRE  SYSTEM 


In  a  tracer-directed  fire  mode,  the  gunner  perceives  both  the  tracking 
error  and  the  tracer  error  on  a  two-dimensional  visual  display.  The  tracking 
error  is  the  difference  between  the  target  angle  0^  and  the  barrel 
pointing  angle  0fi.  It  is  also  referred  to  as  "lag  angle"  later  in  this 
report.  The  tracer  error  e^  is  the  difference  between  the  target  angle  0^ 
and  the  projectile  ending  angle  0  .  In  the  simulation  experiment,  the  pro¬ 
jectile  flight  path  ended  at  the  range  of  the  target.  Each  tracer  round 
disappeared  at  this  point,  0^,  from  the  display.  The  detailed  description 
of  the  configuration  is  described  in  Wei  (1981).  We  will  briefly  summarize 
the  underlying  dynamic  system  in  this  report.  Figure  1  is  the  block  diagram 
of  an  AAA  tracer-directed  fire  system. 


i - 1 

I  blanking  I 


Figure  1.  Block  Diagram  of  an  AAA  Tracer-Directed  Fire  System 


At  any  given  time,  the  target  trajectory  input  is  fed  into  a  visual 
display  device  and  combined  with  the  barrel  pointing  angle  0^,  as  well  as 
the  projectile  ending  angle  0  ,  to  form  error  signals  e^  and  e The  human 
operator  observes  these  error  signals  and  generates  a  control  output  u  via  a 
controller,  or  H-grip,  displacement.  The  control  signal  then  drives  the 
barrel  and  rate  control  plant  for  a  new  barrel  pointing  angle  0^.  Tracer 
round  is  fired  at  this  angle  and  passes  through  the  projectile  ballistics 
computation  to  obtain  the  projectile  ending  angle  0^.  The  task  of  the 
gunner  is  to  constantly  align  the  projectile  ending  angle  to  the  target 
angle,  i.e.  to  minimize  the  tracer  error  e^.  The  dynamics  for  the  elevation 
and  the  azimuth  firing  system  are  very  similar.  In  addition,  the  elevation 
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system  can  be  decoupled  from  the  azimuth  system.  However,  the  azimuth 
system  cannot  be  separated  from  the  elevation  system  due  to  a  coupling 
factor  cos  (eg)£L  in  the  measurement  equation. 

By  introducing  a  state  vector  x±(t)  =  [x^t),  xi2(t),  xi3(t)]r,  "T"  means 

"the  transpose  of,"  with  x^t)  =  0iT(t)  -  eiB(t),  xi2^t)  *  eiT(t)  “  6iP^ 

and  *  0^T(t),  i  =  1,  2*,  the  following  system  and  measurement  equations 

which  represent  the  underlying  tracer-directed  fire  system  can  be  derived, 
see  Wei  (1981). 

x±  *  ui(t)  +  E^t)  ui(t-x)  +  F±  0iT(t)  +  G±(t)  (1) 

and 

Zi(t)  =  ^(t)  x±(t)  i  =  1,2  (2) 

where 


*If  not  otherwise  specified,  the  first  subscript  index  i  represents  the 
elevation  (i  *  1)  or  azimuth  axis  (i  -  2),  while  the  second  index  repre¬ 
sents  the  i-th  element  or  row  of  a  matrix. 


c2  »  cos  01B(t) 

e1(t)  *  -1.34  x  (l-x)  x  [1+  (0.0052t+0.000486t2)  sin  61B<t-T)] 
e2(t)  =  -1.28  x  (1-t) 

gx(t)  *  (0. 0052+0. 000972t)  *  t  *  cos  0iB(t_T) 

g2(t)  =  0 


'T*  ui*  ^il*  anc*  ^12  ^ enote  elevation  or  azimuth  components  of  the 

target  acceleration,  the  gunner's  control  output  and  the  observed  tracking 
error  (lag  angle)  and  tracer  error,  respectively. 


If  we  introduce  a  transformation  on  the  states  x^  and  by  =  cixil> 
*  c^  x£2’  xi3  3  xi3  then  Equations  (l)-(2)  can  be  rewritten  as  follows. 


A[(t)  x^  +  B^(t)  u±(t)  +  E^(t)  \i±  (t-i) 


+  F, 


eiT(t) 


+  G!(t) 


X,(t)  **  C’  x^(t) 


< 


where 


and  x|(t)  =  x^(t),  x^3 (t>J T-  Notice  that  the  coupling  factor 

cos  is  removed  from  the  measurement  equation  and  absorbed  into  the 
system  equation.  Equations  (3) -(4)  represent  a  nonhomo geneous  linear 
time-varying  system  with  a  time-varying  delay  in  the  control. 


Section  III 

AN  AAA  GUNNER  BLANKING  MODEL 


In  Wei  (1981),  the  author  proposed  an  observer  gunner  model  for  gunner 
performance  in  a  tracer-directed  fire  system.  The  function  of  the  gunner 
can  be  decomposed  into  two  parts  to  be  modeled.  In  the  first  part,  the 
gunner  observes  continuous  signals  and  makes  an  estimate  of  system  states 
based  on  his  internal  model  of  target  motion.  In  the  second  part,  the 
gunner  utilizes  the  observed  and  estimated  states  to  form  and  exercise  a 
control  action  in  order  to  achieve  his  objective.  The  former  one  corre¬ 
sponds  to  an  estimation  process,  while  the  latter  corresponds  to  a  control 
process.  The  reduced-order  observer  is  used  in  conjunction  with  a  linear 
feedback  control  law  to  model  the  gunner’s  function.  The  structure  of  the 
model  is  shown  in  Figure  2. 


Figure  2.  Block  Diagram  of  an  AAA  Gunner  Model 

This  model  structure  is  retained  for  the  blanking  case  except  the  ballistic 
equation  is  no  longer  parameterized  by  a  linear  equation  relating  0p  and  0R. 
Instead,  a  more  realistic  nonlinear  ballistic  equation  is  used  as  shown  in 
Equations  (3)-(4).  In  addition,  time-varying  gains  are  used  to  model  the 
effect  of  observation  interruption.  We  will  discuss  the  no-blanking  case 
first,  then  the  blanking  case. 


NO  BLANKING 


The  equation  representing  the  gunner's  internal  model  of  the  tracking  and 
firing  system  can  be  written  as 

i^(t)  -  A^(t)x^(t)  +  B^(t)u±(t)  +  E|(t)  ui(t-x)  +  G’ ( t)  (5) 

^(t)  -  C'±  x|(t)  i  *  1,2  (6) 

Since  both  x^  and  are  measurable,  the  only  state  that  needs  to  be 
estimated  is  x^.  The  state  reconstructor  equation  for  can  be  derived 
by  applying  the  reduced-order  observer  theory  (Luenberger,  1971) 


A  A  •  • 

*i3(t)  *  -(kil+ki2)cix13(t)  +  kilyil(t)  +  ki2y12(t) 

-bikilciui(t)  "  b1ki2ciei(t)ui(t-T)  -  ki2ci6i^t) 

‘  kilCiCi"lyil(t)  ‘  kilCiCi"lyi2(t)  (7) 

The  objective  of  the  gunner  is  to  minimize  the  tracer  error  so  that  a  maximum 
probability  of  hit  could  result.  In  other  words,  the  gunner’s  response  in 
the  control  process  would  be  to  stabilize  the  underlying  system,  especially 
the  tracer  error  therefore,  a  linear  feedback  control  law  of  the 

following  form  is  designed  to  achieve  this  objective. 

u±(t)  ■  (O  +  vA(t)  (8) 

/ 

where 


[y 


il*  Ti2 


»  y131 


is  a  vector  of  controller  gains  to  be  identified, 
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is  a  vector  of  measurable  states  and  estimated  state,  v^(t)  is  a  remnant 
noise  function  assumed  to  be  Gaussian  with  zero  mean  and  a  covariance 
function 


E  v,  (t)  v.  (s) 


i  •  [■« 


+  a 


±2 


« ii  <e> 


+  a 


13 


eiT(t) 


] 


<S(t-8) 


(9) 


for  all  t  and  s.  are  nonnegative  model  parameters  to  be  determined. 

0iT  and  ©iT  are  estimated  target  angle  rate  and  acceleration,  respectively. 
Equations  (7)  and  (8)  represent  the  gunner's  response  in  the  estimation  and 
control  process  of  the  tracking  and  firing  task.  If  we  define  a  new  state 
vector 


^(t)  -  xuW*  *13<t)  -  x13<t>]T 

then  the  state  equation  of  the  closed-loop  system  is  obtained  by  combining 
Equations  (7)  and  (8)  with  Equations  (3)  and  (4)  of  the  actual  tracking  and 
firing  system. 


A^OX^t)  +  D^t)  (t-r)  +  F^U)  +  (t)Vi(t) 

+  Eil(t)vi(t-x)  +  R^(t)  (10) 


where 


*> 


There  are  seven  model  parameters  in  total,  i.e.,  ki,  y^,  yi2,  y^  a^,  a±2f 
and  that  need  to  be  determined  from  the  empirical  no-blanking  tracking 
data. 

BLANKING 

The  optical  display  of  target  was  blanked  periodically  according  to  the 
duty  cycles  and  durations  listed  in  Table  1. 


TABLE  1.  BLANKING  CONDITIONS 


Condition 

Duty  Cycle 
(%) 

Blanking  duration 
(sec) 

1 

25 

1.5 

2 

25 

3.0 

3 

25 

6.0 

4 

50 

1.5 

5 

50 

3.0 

6 

50 

6.0 

7 

75  * 

1.5 

8 

75 

3.0 

9 

75 

6.0 

10 


100 


1.5 


The  duty  cycle  is  defined  as  the  ratio  of  the  blanking  duration  to  the  cycle 
time.  The  blanking  duration  is  the  length  of  time  that  the  target  is 
blanked  so  that  the  subject  cannot  see  the  target.  The  blanking  always 
occurs  at  the  last  portion  of  a  cycle  and  may  reoccur  periodically  over  the 
entire  TOW  firing  period.  An  example  of  blanking  sequence  is  given  in 
Figure  3. 

BLANKING 

PERIODS 


4.0  8.0  12.0  16.0  20.0  T(SFC) 

Figure  3.  Sequence  of  Blanking  for  Condition  4 

m 

The  gunner’s  performance  deteriorates  considerably  under  observation 
interruption  via  blanking  the  target.  In  Yu  (1981),  the  effect  of  blanking 
on  the  gunner’s  tracking  performance  was  modeled  successfully  by  degrading 
the  gunner’s  estimation  gain  k(t)  and  the  controller  gain  y(t).  A  similar 
approach  is  adopted  here  to  model  the  effect  of  blanking  In  a  more  complex 
firing  task.  More  specifically,  the  observer  gain  k^  and  controller  gains 
and  y ^  which  pertain  to  the  observed  states  and  x^  are  assumed  to 
decrease  exponentially  as  the  blanking  starts  and  to  increase  exponentially 
as  the  blanking  stops  (see  Figure  4) . 


Blanking  period 


J  Recovery  j 


(repeated) 


Figure  4.  Degradation  of  Model  Parameters  in  Blanking  and  Recovery  Period 


Given  a  blanking  period  [t^,  t^]  followed  by  a  recovery  period  [t^,  t^]> 
the  degradation  of  gains  can  be  expressed  by  the  following  equations. 


During  the  blanking  period: 


^(t)  =  ki(to^  exP  r  77“ 

T  ik 


Yn(t)  "  Yii(t0)  exp  -  T. 


Y12(t) 


an(t) 


Yi2(t0)  exp  -- 


0^,(0  fl  -  exp  (-  - — 


During  the  recovery  period:* 


kA(t) 


aU(t) 


+ 


°ll(tl) 


(15) 


(16) 


The  time  constants  associated  with  each  gain  parameter,  an<* 
0til^tl^  are  determined  from  the  empirical  tracking  data  collected  in  the 
blanking  experiments,  as  shown  in  the  next  section. 


*In  the  simulation  program,  the  length  of  the  recovery  period  is  defined  as 
the  minimum  of  1.5  sec  and  one-third  of  blanking  period  to  avoid  covariance 
being  negative. 
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Section  IV 

PARAMETER  IDENTIFICATION  AND  SIMULATION 


The  least-squares  identification  program  developed  in  Wei  (1981)  was 
modified  to  identify  the  no— blanking  parameters.  Equation  (10)  can  be  first 
decoupled  and  then  approximated,  via  the  Average  Approximation  Method,  by 
the  following  ordinary  differential  equation  (Banks  and  Burns,  1978). 

W^t)  -  N^t)  W^t)  +  ^(t) 

xi3(t)  =  eiT(t) 

xi4(t)  "  ~k^(t)c^  xi4(t)  + 

where 


(17) 

(18) 
(19) 


+  b1c1Yil(t) 

0 

1 

T 

0 


(i  +  biYi3)  c±  xi3(t)  -  blClYl3xi4(t)  +  btc±vt{t) 

CiXi3(t)  +  Ciei(t)  {Yi3Xi3(t-T)  "  Yi3xi4(t_x)  +  Vi(t-T>}  +  ci81( 


IDENTIFICATION  OF  MODEL  PARAMETERS 


The  equation  which  governs  the  mean  of  states  is  obtained  by  taking 
expectation  of  Equation  (17): 


W^t) 


N^t)  W^t)  +  Mf  iA(t) 


where 


W^t) 


E  |x^2(t)|,  E  E  |x[2 


(1+biYi3)ciXi3(t)  "  ViYi3xi4 


i.(t) 


:iXi3^  +  Ci®i^Yi3  {Xi3^t-X^  “  Xi4^t-X^}  +  Cigi^ 


The  first  and  second  component  of  represent  the  model  prediction  of 
ensembled  mean  of  tracking  and  tracer  error,  respectively.  On  the  other 
hand,  the  covariance  matrix  P^(t)  satisfies  the  following  equation: 

P^t)  -  N^t)  P^t)  +  1^(0  N^t)  +  L^t)  (^(t)  L^t)  (21) 


where 


P,(t)  -  E  )  W4(t)  - 


Wi(t)j  [wl(t)  -  wt(t)j' 


'-'I  •  '  «  *  . 


0 


V*) 


bici 


0 


0 

0 


ciei(t) 

0 

0 


^(t)  - 


- 

* 

ail(t)  +  ai2 

®lT(t) 

+  ai3 

eiT(t) 

- 

- 

0 

ail(t-T)  +  ai2 

e1T(t-T) 

+  ai3 

5 

The  first  and  second  diagonal  element  and  pi22(t)  of  represent 

the  square  of  the  model  prediction  of  standard  deviation  of  tracking  and 
tracer  error,  respectively.  Notice  that  time-varying  parameters  are  assumed 
for  Yil#  y12>  and  ail  to  reflect  the  effect  of  blanking.  Since  the 
blanking  effect  to  the  estimation  of  target  velocity  xi3(t)  is  predominantly 
expressed  through  degradation  of  k^,  there  is  no  need  to  consider  a  time- 
varying  y13*  ai2,  and  a±3. 

The  steady-state  value  of  the  parameters  are  first  identified  via  a  least- 
squares  curve-fitting  identification  program.  The  reference  curves  to  be 
fitted  are  obtained  from  empirical  tracking  and  tracer  data  collected  in  the 
manned  simulation  experiments  without  observation  interruption.  These 
experiments  were  conducted  on  an  AAA  simulator  at  the  Air  Force  Aerospace 
Medical  Research  Laboratory.  Three  simulated  helicopter  trajectories 
ranged  1500  M,  2000  M,  and  2500  M  from  the  AAA  system  were  used  as  target 
trajectories.  Figure  5  shows  some  characteristics  for  the  1500  M  trajec¬ 
tory.  Let  x^(t)  and  x^2(t)  be  the  empirical  ensemble  means  of  tracking 
and  tracer  errors  and  s^(t)  and  s^2(t)  the  corresPond*n8  standard 
deviations.  These  empirical  means  and  standard  deviations  were  obtained  by 
averaging  and  computing  the  variance  of  the  empirical  data  from  40  simula¬ 
tion  runs  with  the  same  target  trajectory  and  the  same  subject. 


dMliiifWtt 
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The  parameters  were  identified  by  minimizing  the  cost  function 


<*.)] 

defined  as  follows: 


i  =  1,2 


where  tg  is  the  initial  time  when  a  selected  tracer  round  reaches  the  range 
of  the  target,  t^  is  the  time  when  the  last  tracer  round  is  fired,  is  a 
positive  weighting  factor  chosen  to  be  one  in  the  identification  runs. 

The  direct  search  algorithm  developed  in  Wei  (1981)  was  modified  to  identify 
the  steady-state  values  of  the  parameters.  The  tracking  and  tracer  data  of 
the  helicopter  trajectory,  ranged  1500  M  from  the  AAA  simulator,  without 
blanking,  were  used  to  obtain  the  following  steady-state  parameter  values 
shown  in  Table  2. 

The  time  constants  associated  with  the  parameters  were  determined  empirically 
from  the  data  of  the  1500  M  trajectory  with  blanking  condition  5  and  listed 
in  Table  3. 


TABLE  2*  STEADY-STATE  PARAMETER  VALUES 


♦Insensitive  for  elevation  case,  no  degradation  is  necessary. 


Notice  that  the  time  constants  for  t,  ,  r  ,  and  x  are  the  same  for  both 

k  y2  al 

elevation  and  azimuth  gunner  model.  This  is  as  expected  because  the  gunner 
manipulates  the  H-grip  for  elevation  and  azimuth  tracking  indlscrlmlnantly 
with  respect  to  the  observation  interruption.  On  the  other  hand, 
and  a^(t^)  f°r  the  azimuth  case  are  considerably  greater  than  that  for  the 
elevation  case.  This  reflects  the  steeper  Increase  of  uncertainty  to  the 
target's  position  along  the  azimuth  axis,  because  the  azimuth  component  of 
target  acceleration  is  much  higher  than  the  elevation  component. 

SIMULATION  RESULTS 


The  gunner  model  was  implemented  on  a  CDC  CYBER  175  computer  to  simulate  the 
man- in- the- loop  AAA  tracking  and  firing  task.  For  the  convenience  of  numeri¬ 
cal  computation.  Equations  (20)  and  (21)  are  discretized  into  the  following 
form: 


w  0+1 


*n  jj11  4.  nn  rn 

-i  -i  — i 


(23) 


R"  Q?  R?  T 


(2A) 


where 


r  - 

r 


■:  ■/' 


exp [N. (t  )  •  o]  do  •  tT 

l  n  — ■ i 


■:  ■( 


exp[N1(tn)  •  a]  do  •  L1(tn) 


2”  -W 

A  =0.06  seconds 


A  simulation  program  was  developed  which  uses  the  recursive  Equations  (23) 
and  (24)  to  simulate  a  closed-loop  AAA  tracking  and  firing  task.  Inputs  to 
the  simulation  program  are  the  time  history  of  range  and  acceleration  of  the 
target  aircraft,  the  initial  angular  position  and  velocity  of  the  target, 
the  number  of  blanking  intervals,  and  the  blanking  intervals  in  chronological 
order.  Outputs  of  the  simulation  program  are  model  predicted  mean  tracking 
error  and  its  standard  deviation. 

Simulation  results  are  shown  in  Figures  6  through  17  for  the  blanking 
conditions  3,  4,  3,  6,  7,  and  10.  The  solid  curves  in  these  figures  are  the 
empirical  data  which  are  obtained  by  averaging  the  results  of  40  experi¬ 
mental  runs.  The  dashed  curve  is  the  model  prediction  of  ensembled  mean 
and  standard  deviation. 


Figure  6  and  Figure  7  show  the  comparison  of  model  versus  empirical 
elevation  mean  and  standard  deviation,  azimuth  mean  and  standard  deviation 
of  both  tracking  errors  (lag) ,  and  tracer  errors  for  the  no-blanking 


case.  Figure  10  and  Figure  11  show  the  results  for  blanking  condition  5 
which  had  a  50  percent,  3.0  second  blanking  occur  during  [11.01,  14.01] 
seconds. 

Of  particular  Interest  is  the  comparison  of  the  empirical  standard  deviation 
in  Figure  6  and  Figure  10.  The  effect  of  blanking  the  target  to  gunner’s 
performance  is  clearly  demonstrated  by  the  sharp  Increase  of  the  standard 
deviation  of  tracking  errors  during  the  blanking  period  [11.01,  14.01] 
seconds.  This  effect  is  very  well  modeled  by  degrading  selected  model 
parameters  as  Indicated  in  the  model  prediction  curve  in  Figure  8.  Similar 
agreements  between  the  empirical  data  and  the  model  prediction  can  be  found 
in  Figures  8,  9,  and  12  through  17. 

These  figures  show  that  the  designed  gunner  model  can  provide  consistent 
prediction  of  the  gunner’s  empirical  tracking  data  as  well  as  the  tracer 
error  data  for  both  no-blanking  and  blanking  cases.  These  figures  also 
demonstrate  that,  for  a  given  AAA  weapon  system,  the  same  set  of  parameter 
values  and  associated  time  constants  can  be  used  to  predict  the  human 
tracking  and  tracer  errors  for  all  simulated  blanking  conditions. 

However,  due  to  the  fact  that  the  helicopter  trajectory  has  very  low 
elevation  axis  maneuvering,  these  parameters  may  only  hold  for  similar 
types  of  low  maneuvering  helicopter  trajectories.  Reidentification  of 
these  parameters  may  be  needed  for  other  high  maneuvering  trajectories. 

The  computer  execution  time  of  the  overall  simulation  for  an  18  second 
helicopter  trajectory  takes  about  5.60  cp  seconds  on  a  CDC  CYBER  175 
computer . 


Elevation — No  Blanking 


Figure  7a,  Mean  and  Standard  Deviation  of  Tracking  Error- 
Azimuth — No  Blanking 
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AZIMUTH  LAG  AZIMUT 
SUBJECT  33  SU8JEC 
TRAJECTORY:  HELIC0P1  TRAJEC 
CASE  4  CASE 


Azimuth- 


Figure  9b.  Mean  and  Standard  Deviation  of  Tracer  Error 
Azimuth — 1.5  Seconds,  50  Percent  Blanking 


Figure  11a.  Mean  and  Standard  Deviation  of  Tracking  Error 
Azimuth — 3.0  Seconds,  50  Percent  Blanking 


Azimuth — 3*0  Seconds,  50  Percent  Blanking 
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Figure  12a.  Mean  and  Standard  Deviation  of  Tracking  Error 
Elevation— 6.0  Seconds,  50  Percent  Blanking 
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Azimuth- 


AZIMUTH  TRACER  ERROR  AZIMUTH  TRACER  ERROR 

SUBJECT  33  SUBJECT  33 
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Figure  13b.  Mean  and  Standard  Deviation  of  Tracer  Error- 
Azimuth — 6.0  Seconds,  50  Percent  Blanking 


SUBJECT 

TRAJECTORY 


Figure  15b.  Mean  and  Standard  Deviation  of  Tracer  Error 
Azimuth — 1.5  Seconds*  75  Percent  Blanking 


ivation — 1.5  Seconds,  100  Percent  Blanking 


Elevation — 1.5  Seconds.  100  Percent  Blanking 


Figure  17a.  Kean  and  Standard  Deviation  of  Tracking  Error 
Azimuth — 1,5  Seconds,  100  Percent  Blanking 
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Section  5 
CONCLUSION 


This  report  summarizes  the  modeling  of  a  gunner’s  performance  in  a  complex 
AAA  tracking  and  firing  task  under  pseudorandom  observation  interruptions. 

The  highlight  of  the  task  is  that  the  gunner  fires  tracer  rounds  without  the 
aid  of  radar  and  lead  angle  computer.  Furthermore,  the  gunner's  performance 
is  greatly  hindered  by  observation  interruptions  via  blanking  the  target  on 
an  optical  display.  A  blanking  model  is  designed  which  consists  of  a  reduced- 
order  observer,  a  linear  feedback  controller,  and  a  remnant  element. 

The  gunner's  performance  is  parameterized  by  the  controller  and  estimator 
gains,  in  addition  to  the  covariance  coefficient  of  the  remnant.  The  effect 
of  blanking  is  modeled  by  degrading  these  gains  and  coefficients  as  a  func¬ 
tion  of  blanking  duration.  An  exponential  decay  form  is  assumed  for  these 
parameters.  The  associated  time  constants  are  determined  from  empirical 
data  collected  in  the  blanking  experiment.  A  direct  search  method  Is  used 
to  identify  model  parameters  systematically  while  minimizing  the  least- 
squares  error  between  the  model  output  and  the  empirical  data. 

Computer  simulation  of  the  proposed  gunner  model  show  that  the  model 
predictions  are  in  good  agreement  with  empirical  data  for  various  blanking 
patterns  using  a  typical  helicopter  trajectory.  These  results  demonstrate 
that  the  model  can  adequately  describe  the  gunner's  tracking  and  firing 
characteristics  in  an  AAA  weapon  system  subject  to  observation  blanking. 

This  gunner  (blanking)  model  has  been  incorporated  into  the  MTQ  series  of 
POOl/OBS  AAA  engagement  models  and  is  designated  as  program  P001/0BS  3/6B. 

This  composite  program  POOl/OBS  3/6B  can  be  used  in  the  evaluation  of  air¬ 
craft  survivability  and  performing  weapons  effectiveness  studies.  Documen¬ 
tation  of  POOl/OBS  3/6B  is  in  preparation  at  the  Air  Force  Aerospace  Medical 
Research  Laboratory  and  will  be  distributed  separately. 
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APPENDIX  A 

LISTING  OF  PARAMETER  IDENTIFICATION 
PROGRAM  ELEVATION  CASE 


iirfrjaun 


3SIT5 


A1  TACHtlt  P- 1  •  1NU6H  EUSU6  JS  3t  C  Y  •  1  •  NR*  1. 


DIMENSION  ALFHAt7),PSI«7l  ,  A  (7)  ,0(7)  •  ALP  HAW  7) 

0l^Ff»Sf5N_5up5(6riTir}7^<T,7r,r  t7,rj,  XI  (7, 7),  3 (7, 7), F 1 6)  , :?i  6,7) 
CCNMON/ARRAV/XI  1000),  SU000 >  ,EL0D  (lboO)  tTXU0  0  0>  »TSI  1000)  ,  R  ANI 1 00  0 


•  ■  nti  i 

CCMHON/S/CO, CELtNSTP, NOIN,TO, IPT,Y10,Y20 


INTEGER  QfSUEXT »FR 


00  1  1*1,7 


READ',KI,NDI7,TO,IPT,Y100,Y10 

PRINT  4j,xi“,  NOiN,TO,IPTjn® 

FORMAT  T1H1, "NO. OF  PIS  •  “,  14,  2X  .“ORDER*  “,I2, 

-c_«y7**i'!nr  n Ft*  -,si?.5//ix, **Rt«r_Evm_%r?r'r_ooiNr vi o*  mu 

c.si 


DEL«1.03»IPT 


K*K1-XT 
KEAbti,*: 
C  1*1,  XT) 


L.CZNNiTlJUMZ,  AZSO,  01)  M3, 


READ  (1,  ,3)  T  ,bUN6«  AZ,  AZD,  AZOD«EL,ELD,  DUM.AZNN,  0UM2,  AZSO,  0  JM3 


IMtJMUl  4Y* 
FORMAT  ( 12612,5) 


tfwjwjjjTjt 

LmiODj 

mmng 

i:  8  §  cj  m  *  m  1  **  ^ 

HsnnEa 

fiTOiirmd 


T  FET  ( IH)*EL 


IF(F0r <’>)*0«48 


IF(MO0(I-itIFTUN£.0l  GO  TO  47 


_TXtn>«2UH4 

COMTIMUt 


"TTSTWTH^I - 

*E4 ON  (ALPHA  <l)9X*l97) 

"RE»SF»S - 

AJ»0.0 

xnnrmrr 

fc(ii«.i 
"1 0  lrJa07(J 
A(t>*2.3 

"cnrrwua - 

n»so2ir«o  • 

" V  2T"«"CTACT=1 1 nJ-TlTJ 0 re 07TTO  I*T 5 ;FN(*N f  i)>  OrWS* S A WTIT"2 1  •  CO 
i  S(ELO-tlOO) 

“CALI  INTCTXL  FHA  t  AJJ - - 

QLOJ*AJ 

PRINT  4tJ#ALPM - ~ 

„  40  FORMAT  ("OALPHAa  G 12 .5*  **  ING12  •5t**t**96  12.5# "IN  G12.  59-|~i 

*  C  tr^*^7"(,r»GlZ#5>t"I''f  512*  SI 

*  PRINT  41, AJ,  ITfSU6IT 

.  “Cl - nTTKATT-  J*  "»  G12*5t"  irEPATI0H5^7T5;^Wm**TTUNC«-"N 1 51 

<.  11  00  100  K*  19 7 

„ - siFUTCwrrvi - 

CO  12  L*lt7 

*  - Ai»HAMru‘»"ALPhA(u»n  nr^vncm - 

.  1F((L«G£«  21  •  AND*  (L*LE  »4) )  GO  TO  12 

. - IF  (ALPHARlLT7LT".0.in  AUPHITNIU^An^RTa - 

w  12  CONTINUE 

—  CALL  I  HI G  CA  L  FH  A  N TA1> - 

IF  (AJ.GT.OLCJ)  GO  TO  20 
OLJJ^D 

PRINT  40,ALPHAN 

PRINT  41t~A3f  IT «  SUB II - 

U  «l  *  C  (Kl  ♦£(  W 

VT0^3*TPO - 

GC  15  H*l,7 

A  J*HA ( HJ  *TTLP  HAN  (HI - 

C  NTlNUc 

if  imcmrr  nrsr  jmmnro - 

GC  TO  25 

rK^^mn - 

IF  f  A ( K) *L£  »  1*5)  A  (  <)  «0. 0 

xicnrro - 

CG  33  L*l»7 

IFIATU.U.O.MGU  TO  30 - 

CM  1#  0 

CCNTIROT - 

IF  CC*.  NE.0. 0)60  TO  100 

GLNUSU^OtO - 

00  32  M*1  *  7 

'svn^nrri^xi  - 

bLN2*SUN2*XI  C2,lO**2 


lb 


~Zfl~ 

75“ 


,  ^T 


X1*SQ?T(SUH1) 


X21S  JRHSUH2  J 
X3*X1/X2 
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It  X3“ 


FORMAT  (• ljs  •»G12*57»  ALPHA*(**Gl2.5f*lN51?«5tMNG12«5* 
7^l*i3ir,59*l^l2;5fMS;i2,5trlFf0«#5/ -  - 


3+  XX<1>«  NC12.5/*  XI  (1) /XI < 21  *  NG12.5) 


GC  T0~110~ 
100  CONTINUE 
TJT5 - GCTTO”II“ 


f  (L )  a  G  •  3 
F21  l*M)*0.0_ 
CC^Y  I  NUfS 
cc  it?  i«i»7 

00  117'  <«1,  7 

xi< j*k>*o.o 


0C  120  I«l* 

sui3»su^y»(A6S(  o  <  in  i _ 

IF  (SJMT.LE.EFSIGQ  TO  1000 
00  130  *»1*7 


DC  130  J*  1*  7 
ZHt  JlOfNJ*  V<N,JI 


CCUINUE 
OC  140  J* I*  7 
OC  140  L*  1*7 
CC  140  <* J*  7 


XlUtL>*XlU»LJ«ZK»0 
C  t*4TI  NUE 


S114«0.0 
OC  150  J*l*  7 


$UM4*$U14fXl  (1*  J)#  *2 
Sim*SQRT  ISUF4J 


CO  155  J*  1*  7 

V(l*  J»*XI  (1,  J)/SUM4 


K  C'JNf  *  I 
G*<0UNT-1 

otf~  i?§~  <*i7q 
OC  ISO  L*l*7 


F  2(  3*  (Q) #V  f  K*  HI  ♦  F2IQ*rtl 

CC^TtNUi 

F  I3)*_0jiJ _ _ 

CC4ftNU€ 

CC  11G  1*1.7 


nmptni 


20  0  SU15  ( 0 )  *S UH5  <Q)  ♦  8  (K  OUNT.  Ml  **  2 

- sram  *  ^  nsu**T(<3  n - 

OC  215  4*1.7 


IF  IOlHT.LE.7j  GO  TO  159 

-TFrmi - 

SUSIT  *  0 


«• 

OO  250  <*1*7 


1400  CALL  EXIT 

fcM 

SUBROUTINE  I  KT G  (ALPHA  *  EJ) 

, - rrw  otts  mnaroo  jTEoano  ooittx  <10001  ,ts  <100  or 

1,TAU(1030)*TKT<100  0» 


r«vc o< n n ?k# <jn ? tin f  r%s iinrJij 


0I1ENSI0N  H(4>,P(4,4>  ,Pl(4,4),P2(4,4),A».PHMn,A(4,4t 
l,6<  loJ",?BTlTi7EbWt  16T7E  1*741 ';mv,  CQC<  4, <.1,0(41  — 

2,  XI 1 10  9  It  ,X4  <10  00),EOH(10  30)  ,60  0H(  10001 

cwr VALENCE  (4(1,1) ,lFf  I t)V(E Alt ,1 > ,E6 11  >  )  , » £  II  tin  1 » 1)  ,tB  1 1 T  <1 )  ) 
U  IT  A  *T'l./ 


r 


C  INITIALIZATION 


SCAL»C0‘*2/D£l 


OC  1  J* 1 , NO 


CCII)JM0. 


NGI*N0-1 


NAAMD/2-1 


N(I)»0. 


FII>»0. 


P  11, 11*0.0000256279 


X3(  1)  •  0*000 


EMI  11  *X3  111 


IFIARG.JT. -200.1  St*EXPIARGl 


C  COMPUTE  AND  STORE  STATES  X3  AND  X4 


2FI I TOOIXK,  1  00)  »£Q.  0)  .OR.  IKK.EQ  *1 1)  PRINT  9j,T,X3(K<) ,X4IK<1 


9b  FORMAT  14612.  5) 


K2*«-l 


IF|Al*Hlllt  •£<).•. I  GO  TO  4 


►  til  AlOili  LlLHX^'l'lLilUlfQ^  lliilUilAtll 


GO  T0  3 


f  HI  GH  I 


CONTI NU 


C  COMPUTE  ANO  STORE  ESTIMATED  TARGET  VELOCITY  AND  ACCELERATION 


£Mt<K)*X3IKK)-K4(KK> 


N1*N9**2 


C  START  INTEGRATION  LOOP 


M  l»l*T20 

“ISET*fl - 

OC  100  XK*1»  NSTP 


««■  4  m 


IFIKK.LE.H)  GO  TO  100 


IF! ISET.NE. 1 )  GO  TO  101 

sr*  A9i  ihu-xempirki  i  **z * wmsTimi  **  r - 

SS9*ISQRT|P|  1,111-SEMPIll  >•*£♦<  SORT  IP  12 ,2I» -TSI 1  ll*»  2 
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till  Aiut«  ♦  001*1  *«2*TAU(«)  ♦•*»***  )AU| *K)»*  2) »il Ml  THAI  IKK -H>> 

ALl*AlfMA<2> 

- RRWTI - 

Kl*«U 

. - SftNAr/nuRKl - 

A2*1«*(TAU(K1) -TAU(KK) >/oel 

- - CTR-SD^LS - 

A  (1 1 1)  *-C  OR 
4ni7TMT-CTl*JtPHJV3I 
r  A(2,MCl)«A(lti|*AEl*A2 

.  a  i?TNTnnn  rray  *  *a»  a? - - 

„  00  2  1*1,  ND2 

..  *  1*  I  ♦  2 

A IJlt I ) *RA 
A  U1* 

..  2  CCHTINUE 

C7iUrTT3RTTTinfl7^l)IL7tB7tBITn75y - - 

CR1*C0*ALPHA  (41 ♦AL1#A2 
.  C R4»tO* AL^H A  (41 

SCAL1-S3AL*  <AL1*A2)**2 
-  CCTTI,  2 »  *AlPKS  (? f*5CATI 

.  C  IF((HQO(KK,  1 00)  •  EQ«  Q)  *0R«  (KK*€Q«1 ))  PRINT  99,  T  ♦  W  (1)  ♦  Y 1 0*  SHE  ANf 

-."C  i  mriPtUUUSSO  '  —  ~ 

.  9*  FCRMAT  (5G12*  5) 

.T 

.  C  COMPUTE  ME  A  H  TRACKING  ERROR 

•  Z - 

OC  110  I*  1,  N Cl M 
5(TT2tf~" 1*1,  NCI  M 
G(IMGCI)  ♦EA  U*J>*4(J> 

TJ  CCNTINOE  : 

„  110  CCNTIN'JE 

,  FTIT5T  r.^CRn  F X 3TO(K >TCRV* X¥IKK1 - - ~ 

F(2t*K3(KK|  ♦  CR3*(XMKK-Ml  -X3(KK-M))^0*0  01»(  1. -A2>*  (5  •2*0*97  2* 

„  C  TAUURIT^CT^TBFnW-TOl - 

v  OC  130  1*1,  NCIN 

. - cri^ir  j«t,  2 - 

„  C(I)*C  (IMEAINT  II, J )*F ( J) 

.171 — czr rmz - 

«  M(t>*Dm 

*  - - 

«  130  COHTINUE 


„  C  COMPUTE  ERROR  DUE  TO  MEAN  TRACKING  ERROR 
«Z - : - - 

SUSAN* SMEAN*  (M1>-XEMP(KK>)**2MH(2)-TX  <KPU**2 

- 

„  C  COMPUTE  COVARIANCE  MATRIX  _ _ 


CCD<1, 1>*  (ALPHA C5)  ♦ALPHA  (61  MBS (ECH(KK)  >Ml*M  A  ( 7 > MIS  (EOOH  COO >  I 

■  .. - nwi - 

f  „  IF(«.GT.M)  CQC(2,2>*  (AlPHA(5>»ALPMA(6>MB$(EDH(KK-Mn 

f  „ - i  ■ —  A L  p  H  l  <  7  *  *  m  lEOCWK^* n  >  •  SC  At  1 - 

\  GAEL  HULTCEAlNTtCQCtNOlN,  N1,P1,10) 

|  „  CAll  hulT1ET*P*ndih*hi,P2, lffi 

OC  220  1*1, NCIN 

r  . - G«r7?o-j*i,NciH - * 

p<r»ji**icx«  j)»P2<itU> 


iUmjTINE  MUT*E,F,L.Ll.H,rlR) 

TiiwirnjirmT^FUTYsrisi  ,na»~ 

CO  10  I*t,L 


00  11  K’ltl 


OC  5  J»XtUtL _ 

t«p»ttip*e  i  j»*f  an 
ii>n»i 


■  rn  mini 


HlKKt  *T"MP 


IF1MR.E3. U RETURN 

Tjc_?o_r»i7i: - 

00  20  <*ItC 

Tn»ir. - 

1I*K 


Tt<P*TE  1P*6(  J)*E(II) 


SlVDUTINE  DSCRT  ( NQIM *  A«  DEL  »EAt  EA  INT >  NT  ) 

- orrriTsnN  « c  n ,ea(iv, cainthi  ,c  tef?  jot - 

C  SETS  EA*EXP|A*DELI  »E  AX  NT*  INTEGRAL  EA  0  T3  3EL 

- NCTKiiTnmFi - 

NNaNDIN'NOlH 

- nuiawt-i - 

C  CEC  ( NT)  *1, 

- 00"13-ITTTHm - 

II*NT-I 

- r»  CUEFTI  XT  *XJEL*COEFTTI*  1)  /FtOATtI  X - - 

C  NT  MUST  BE  AT  LEAST  3 


00  60  La3,NT 

- CALr-wjtr  (ATE#ii»r,woiH,tiir,  tmi - 

IFfl.EQ.NTlGC  TO  70 

~E0~CALL  DXACnroTHtEAINTrEA»r70tCOEFTL11 
70  OC  10  II»  1«  NKt  NOXHl 


60  CCNTINUi 

- tn - 

blBp3  JT INE  0IA6INDIH,  A,J,C1,C2I 

- ormiioN-*  t  n  ,1111 - 

NClNl*  N3I Ma 1 

- RK*«rarH»NOin - 

NM*NDH*1 

- II*  t - 

IFtCl.EX. 1.0)  GO  TO  10 

- 3C  5  J*I,NN,  K)I*I - 

Kaf'Wl 


MJHM1 
00  6  IaJ.K 

rmnnn - 

*  <111  ■  Mil)  +C2 
m«ir»H3m 
«eru^N 

- ETP5 - 

535,4, 2. 45, 2, -0.017964, -0.020511 

.  1T54J1,  .017491,  .I?443J».423i0,.224£6T~-77iT7 97 5E-J,.  17  302E-~S' . 

.  0.01 
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APPENDIX  B 

LISTING  OF  PARAMETER  IDENTIFICATION 
PROGRAM  AZIMUTH  CASE 
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,  *¥FfftT70Ut  :^7 5  TO  3  V — L760295VWEI  •  256-39  60 - 

.  M«P(QN) . 

,  COMMENT  ;  ’ 0WVFA7 1 C ,  T  0  =L7  B  U  2  9 SVCT  *  1  T^TOEMT- A t '  PA?^^/Sr“HETr~TRA^*  * - 

,  F  l  N# 


.  AT  TACH,TAP£?,TWL6T  RACtRS(J3J  33,37*1,  ID *L76 0 295 ,MR* l . 


PRO 3*AM  OPT  (TAPE l, TAP E2f  INPUT tOLTPUT) 

~  CHENS  I  TN  A L  FHAT 71  VPS  1(71  ,  A  (7T70(7)  ,  A LFHATrm - 

uHEMSI  ON  SUHi,(b|fE<7|  ,*{7,71  ,Z  (7,71  ,XI  (7, 7 1 . 3  (  7 ,7 1 ,  F  (6)  •  •  2  (6,  7  ) 

- tCTH^/A?rA1  /“XVlOOTITVS  (TO  0  0)7*700  (10  DOT,  CISC  rO  3  0)  AN  (100  3) - 

1,  TX(lOai)  •  TS  (1000) 

CCIHON/S/COt  CELtNSTP^  NDINfT0Vn»TVY10tY20t  A7O0 - 

lCGIOAl  c  Mf  P  AR 

INTF3ER  9 «  SU  EIT « FR - 

L  G*  0 


liPWrr 


Nf A  ? 1  *  NJA  R-  1 
- LC  1  I*l,NPAF 


£ C  1  JslfNPAS 

“w  <r,  j  >  -  <  r  (i*  j)«ovo 

IF  (I.ET.  J)  V(I,J)=XI (I,J1*1.0 
~2  irrjr*T,  o - 


bet,  jj* 

i  — ccntiru:  - - - 

F  R*  1 

- R  £5  0  >,  <1,  NOTFtTITfTPTVTlOD'rYl  0 - - - 

PRINT  4?, <1, KDIM.TO ,IPT,T 10 

T5T2 - FaMXTTlHI#-Rt!;or  PTS~*  » ,  I  v,  ZXV^OROEfrr-^Tt? ? - 

C  2X  ,  ”  I  N I T  TI7E=  ”, Gl2,5// IX,  ”RE  40  EVER  Y  M?rH  p  OINT  Yl  0*  ~,S12 

C.51  .  -  -  ~  - 

K7«T0Z  3.03 

K  *<  1  *<T  - - - - - -  '  ‘ 

OtL*3.31»IPT 

- TTTT  =  C - 

nEAD(lt  *3) (I ,CUM4  ,  A  It AZD, DUM1, EL, ELD, EL CO,D JM2, ELMN,  0UM3t* LSD, 

- C  1*1  V<T)  '  ‘  - - -  -  -  -  •  -  •  •  •  -  •  * 


uc  >4  ih,k 

- REACH,  *3)  7,CUHL,  AZ,  AZD,-OUMl>Tt,£tO,ELTJDVOJM2,rLMN,  DUM3,EL3D 

lFCEJ-Cin  *,4,45 

TT3 - FORMAT  (I2G12.5) - - 

Jo  lFC*n:(I-l,IFT).NE,0)  GO  TO  44 

Ih*  (I-D/IPT  *1  - - - 

1F(  I-UNI*  t)  GO  TO  49 

AZ10*6Z) 

1f9 - CONTINUE - 

A  IH»«DUM1 

X  (IH)  O  JM2  - -  - 

a>(  imuEl-elnn 

■  . - -  S  (THI  O  JN3  - -  - — - - 

n  A’M  t  HI  *7 ,5 

- IFrDUNJ,,LT;2877V-T-RAN(IH)xOUt14r(930T^Tl"9M)Ji4) - 

M«?fiN(H)/jEL 

- *  1 F (  H.LE.  H)  GO  TO  44  - — - - 


i  T  a  I  SET ♦  ! 


-  IFfISET.Nfc.1  )  GC  TO  44  “  - 

a  ZT  a  i  *  A ! 

vi - conti  n  ur - 

R  fci  ?  (  ? ,  -*6  )  (T  lM£,aUN5,  TEL,  0UM7  ,T  ELSD,  1  0«  ltKH 

LC  4T  I  »1 ,  K  . . . 

K€>: (?*  ♦&>  T IME, GUrt5, TEL,  OUM7,T  ELS3 

*  ~  *  IF(P?r  (?)  )  46,46  - - 

4b  FC*wW(5Gl2.5) 
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4* 


r 


4«  1FM03II>1,IFT).NE.S.I  SO  TO  47 

IG*U-1»/IPT«1 
„  T5fTT5r»5UH5 

..  TS<IG)»)UN7 

.-57 - ccNTisur - 

.  CO*  1.7# 

. - Rsr  p*r>^T - 

htAO»,  ULPHA  <I>.I»1.NPM> 

.  S*S 

.  POINT', EPS,  (ACPMUI,  1*1,  NPAR) 

.  *  j*0.0  -  ----  “  "  "  "  . 

„  10  OC  8  I*l,NPAfi 

„  t<IT».l 

0  111*0.1 

-  A(I>i?.U 

.  •  CONTINUE 

.  n«swzr«v 

V2<M*ZT1U>UZ0>V100I 

. - LILT  I  ^fG  iiL  PhTTJ  JT - 

.  CIS  J*A  J 

-  PRINT  4 f.TSLTHTrNPrTNP^l.  NPAR) 

..40  FC2NAT  (-0ALPHA*  “, G 12 .5, " l“,  G 12.5,  "l“,  6  12 .5.  “ I ",  G12.  5.  ■ 

»  6  C12.5*“I“,G  12.S»**I  “»  G12.  51 

PRTNT  41,  AJ,IT,$U6IT 

.  -51 - raf«rfi»'T«  ?7 5  r^ntT rsATiows*  -7r5mninorrTroTsir 

.  11  00  100  <*  It  NPAR 

.  JCSZT*S'JSITT| 

„  OC  12  L*1 , NPAR 

. - AitfMH.ra*“iL  fH«n  »n  kj  *7Tff,n - 

«  IF((L.GE.?I.AN0.(L.LE.4)>  GO  TO  12 

« - irTTTwrni mm  .a roi  'alpha wrns-ALTHunn - 

.  12  CCNTINUi 

.  cwnriNTG  (xtwANTJjr  ~ 

.  IF  »AJ.  JT.OLCJ)  GO  TO  20 

. - clu*»i - 

«  PRINT  43, (ALPhANINP) ,HP*1,NPAR> 

.  PRINT'  41,  AJ,  Ii, suSTf 

C«)*3(<)  *tiK) 

«  EW»5MC»CI 

.  OC  15  S*l,NPAfc 

„ - XIPflSTHlVAlP  Bnmo - 

.  1*  CCNTINUE 

. - inr«i.cr  ;i:?r cKirrm - 

..  G C  TO  25 


rMrmtlii 


IF  (4nO.LE.lo9)  A00»0.0 

“CMffTd - 

CtC  31  L*1 .  NP  AR 


CMliO 


105 

110 


bC  TO  U 
SHJ=0.0 


~C  C"  f  l5'C«17KMar 


GC  115  N*  1 ,  N FAR 


TCUTiOTT 


.“rrr 


F2tL#K)=0.0 


Trc^rrrjr 

GO  117  J*  It  NFAR 
"DOiir“<»  1,‘NFjnr 


.  117 


XIC J,K)=*0.0 

~ccT?frr«itNTs?r 


120  SU3=5U13MAeS(0(I>>> 


“i  f  csTrnrcETtrsTGD  ro  nnnr 

GC  110  N*  1 1  NPAR 


TC'nO~~J«  l*>iFAR~ 


Z  <Nf  J)  =0(N>*  V(N*  J) 

cnttnue - 


GO  140  i=  1 ,  NFAR 


lbO 


"DC  1LD  LkliNFffR - 

OC  140  <* J »  N FAR 

"xt xi  CT,D>zTxvrr 


CCNTINU! 

"SLH44(T.a~ 


LC  150  J«  It  NfAR 


.  "mr 


"sut  sm  c  >Tr r  rr J  r*  *■? 

SUbaSCW  (SUf-J 
"TCT1 


~J«T#  N  FAR- 


1*>5>  v  (l,  Jf  sXKlt  J1/SUM4 


159 


“RCJNT  *2 - 

b*<0’JNT-l 

TJcrrrr-OTVir 


lbO 


GC  160  L  =  It  NFAR 

7m«fn)  ♦  xrtKouTrrrcr^rncrrr 

cqntinje 

- DC  170  HbITWFAR - 

F2(3,  M)*F(Q)*v<KtM>  ♦  F2(dtM) 
T7TJ - CCTTINin - 


FO)  =  Q.3 

'continue 


*  T9TT 


GC  190  1=1, NFAR 

"B  K  O  J  NT  f  D  «  XTTX  CUNTVTT-F2  COTIT 
SU15  ( Q) =0  •  0 


itnfar" 


2u  0 


SU'lptD)  *SUM5  <Q)*B<KOUNTfR>**2 

“SDi5ta>  ^sort  (suHsnn - 

GO  215  1=1, NFAR 


- V  (OUNT ,M)  *G  (KOJNT,  NT /SUN5TQT" 

215  CCNTINU" 


"KCTNT  fOtJKTfT" 


XF  (OUU.Lc.NFAR)  GO  TO  159 

~rr*rr*i - 


bU3It=3 
Ffi*  1 - 


CC  250  <* 1 , N FA R 


TKirrr 


OKMC.3 

“iiKwrr 


250 


cunt  i  no- 

- GC  TO' it  *“ 

IbOO  CALL  It  IT 


:n 


SL15)UTINc  lKTGULFHAtEJ) 

C  CNMON/  IF  FAT  /XcMP  ( 1  00  01  t  SEMPTHTOOlirEOO  (TOOO  >  VttG  (100 Dl  ,  TAG  (1  000  > 
1,  TX ( 100 01  ,  TS  (1000) 

-CClFON/;/CO,CcL,NSTP7NO,TariPTiTlOVY20r*2OO' 


CHENS  I  IN  RU)  ,P  <4t*1  ,  Pl(  4,  *)  fP2(4y  4)  ,ALPHA(M  ,  A  (btt) 
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1,  (lib)  •:B(16).£61NT  U6),:Mh,<,)  fEAl  NT  <4 , 4>»  -  U)  ,  CQC(  4, 4) ,  C  <4> 
2t%U  1 0^3  J)  «X4  <10  00)  t  EON (10  001  .EOOHUOOOI 

~ttfJI*iC*NC£“  (A  ( iVil  t 9  (  i)  l  t  <t  A  (l  v  1) «  £B  UlFi'dAlHr  <1*1 J» EBlNT  (1)1 
OAT  A  WT/1./ 


*  C  INITIALIZATION 


c.e*co5<;lg<kioj 

Afc js-GEC*  AL  PHA"{  ii*C 9 

IF  C  A  ^3  *GT • -  2  00* I  S 1»EXP( A  RG) 

X  3(  K 1 ) *X3  <KK)*tDD<<<) ♦OCL 
X F { AL?HA( 1 ) •EG*  0 •)  GO  TO  4 

“x  4  <  mss  i*  nr  40?k  j  ♦roorm  •  ran 

GC  TO  1 

X  4fK  l f »  X  4(  KiTReCO  f*  K  )*  o  £  L 
CONTINUE 

fc  fa  <  i  t  ^  *ti  tnr-x  v  r<  1  > 

1  F( << •  GE«  2)  ECjH  CKKI  «(£OHCKK)-EOH<KK-1M  /OEL 

"7i*T  ATJ  (XO/OEC - 

IFC«.LE.MI  GO  TO  10 
~I$ET*rS£Tn 
IFI I3£T«N£»  1  )  GO  TO  10 

“wm*ri3*ca 

h  (2)*r23»C6 


KHUM 


P(?t  ?>*<0  •00  71567*CB) *•? 

"s>-  t  ^  ru i  J  /  ca-  xenp  rm  ^2>~cw(E>/cB-rxTTn  ”  z 

SSO«(5^r  (P<l»in/CB-SEHP«1))^  2HS<m<PC2»  2)  I /CB-T3 1 1)1** 2 
ccitinue 

“7*70 - 

DC  100  <K  =  IS Tf  NSTP 

kI*«*x 

M«TAJ<«)/OEL 
RA*N AA/TAU  <  K Kl 
OC  2  I*  l»  NO  2 


T6*-T*lEiO*TA  M£l_G(KK)  I 
SCAl«(C3*Ctil**2/0EL 

~T«rott - 

A2»l.-<rAmKl)-TAU<*K))/3gL 
”CTR*C0MLPHA  «2)*C8 
A(ltl)»-C0RMb 


Ad#N0l)«-C0l»*A2 

- rn7Nrj*nit2)*»g - 

A (?*  2) »T9 

X 

C  COMPUTE  TSAJSITICN  MATRIX  £A  AND  ITS  INTEGRAL  EAINT 

"C - 

CALL  0S0RT(NCIHvb«9ELtEB«E8INT95) 

t4V*CB - 

C 

X“ START'  INTPGRATIOITUJOP - 

C 

X - intHDOIKK,  10J.Ea.1i;UR.IXX.tU.  1>>  PRINT  W*TfW1 

c  i  stptipii, n  i«sso 

"99 - FTTMATTSGT2.~5) - 

C 


00  1?0  J«1»NCIN 

- cnmrm  *ti(rnrTm - 

120  ccntinus 

tto — cc^riwui - - 

f  (IMIC3-CR4)*X300C)^CR4»X4(KK) 

- rtn^xinocr^XB - 

IFKK.GT.MI  FI  21 »F (  2)  ♦CR4i  <X4(KK-M)-X3(KK-M))  M? 

- ocrutr  i»irKnN - - - - 

OC  140  J« It  2 

- DfII*DTI>  ♦  EAINT  lit  JI^FTXl - 

140  CONTINUE 

- n*i n  *  cm - — - - 

OCtMO. 

TTO - CONTINUE - - - 

c 


- yprmyiEHN»  wm/CB-txgwrtKTrr»»2»m2ivcB»rx(Kr»rsm»»2 - 

c 

C"C0HPUTETTCV1PTANT*TTUTRIX - 

c 

- C C3  I  ITI)  «  fACFHri5T»AL1»lfArGr^A0STrtJHt  KKn  ♦ALPHA  )*A3$  CEOOH  ( KK)  I ) 

1  •SCAL 

- CCTT2t  21  »ALPPAI5)»SCAL»A2» ♦? - - 

1F<<<*GT,M)  COC (2* 2 )* (ALPHA IS) ♦ALPHA I  6) * ABS I  EDHf KK-M) ) ♦ALP  HA  17)  ♦ 

- 1  3SfEO:)HfKJC-Mm*$CALM2**r -  — 

CALL  MULT  (EA  INT  #CQC  #NOIM»  Nl»Plt  10) 

- CALtr-rULTCEAfP#KOlHtNlt>?»10) -  *  -  * 

OC  ??C  I*1,NCIN 

- CC'"?20  J»  1 1  ACXH - 

Pitt  i)«Piu«  j)^P2il*J) 

•m —  cc*ittnue  '  - - 


)-TSm-IST))“ 


^COMPUTE  ERR3P  OUE'Tff  STANDARD  DEVIATION 
C 

- 537*>50VCSQRT(P  fit  1)1  /Cr^S cnrOTUT  +* 

A**? 

nrr-  continue- . - - - 

£3*  0£l  ♦(SHEA  MtN  T*S$0) 

- RETURN  . . . . 

EM 
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V AfU tLt tflXNl  «  Nl  I 

0  Ii£NS ION  A  <1)  *  EACH  •  GAIN  TUI  #C  CEF(30) 

T“  S£TS  £A*tXMA*OEL)  t£ArNT*lNTcGRAL  EA  ITTy  OEt 

NCt  *1**01  *♦! 

NN*NOlH*NQlH 
Nl  'll*  NT  - 1 
CC£FTNTf*l. 

CC  13  I*ltNTPl 
Il*NT-t 

10  CC£-  (II>=OEL*CO£F<I  1*1) /FLOAT  U> 

- - NT“MUST~6£~  ir~LEA3r3 

CALL  DUG  (NDlM*tAlNT,  A  f  Z0  EF  ( 1  )t  COEF  ( 2 II 


OC  b 0  L*3?N1 

CALL  KUlT  ( A  *  £AI  Hi#  NOI  M.NN  *EA  »  1) 

- 1  r<  l  t  *  cr  t  it  7  0 

bO  CA.L  CnG(NUlM,EAlNTf  EA,1.0,COEF<U) 

“  n  CO  P  flr  It*  It  N Kf  NOlMl 
EACH )  *£A  (ll)*l«0 

— strcr^rnnT-: - 

ENO 

- SL1  To  uTTti  roTTCTK STOUC  C 17C  21 

0  IiENS  I  3N  A  (  1)  •  B(  1) 

Net  *!•  N3l H*1 
NMNOI  M*NOI  M 

- n  - 

11*1 

i  f  ( cite  iriTff  rwT0~ro 

to  5  J* 1 »  NN t NOXH 
fcC  A  I  *  it  K 

- ir*  irTin  *rm - 

A  (II)  >  A  (ID  *C2 
1  lI*II*N)IHl 

_ RETU*N _ _ 

10  OC  7  J«T,NN,*KOlH 
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APPENDIX  C 

LISTING  OF  AN  AAA  GUNNER 
MODEL  SIMULATION  PROGRAM 


HSM,T20,CM70000 •  L760 ! 95 ,HEI, 25 » J960 
COMMENT. •NEW0WLSTMU6, 1 0*L 760295,  GY* 1* 

CO HHENT.  *AAA  NOOE6  BLANKING  5IMJLATI0N  PROGRAM* _ 

ATtACH,TAPEl,OWL6HELSUBJ33,IO«L7»0295,CY«l,MR*l. 

FT  N. 

L50. 

PROGRAM  SINJSmPJT,  OUTPJT,  TAPED 

COMMON/ S/CD  (21 .DEL •  IH,NOIN,  YiO  (2),X3(2),  ELtELOO,  AZOO.MTAU,  RA,  A2.NO 
A1 ,N02 ,NAA,UEL,UAZ, ELTR.AZf R,ISET.Z1.Z2.TAU.TS(15).TE(15) .T.IBL 


c 

c  THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  SIMULATE  AN  ELEVTN  ANO  AZ  IMUTH_TRAJ:k 
c  TASK  IN  THE  TRA  IE  P-0  I&ECTE0  FIRE  (MODE  6>  SYSTEM 

C  SUBJECT  TO  OPTICAL  3LANKING  _  _ 

C  INPUT!  THE  ELEVATION  (EL)  1  AZIMUTH  (AZ)  ANGULAR  ACCELERATION  OF 

C  TARGET.  ANO  BLANKING  OURATIONS(UP  TO  IS)  IN  CHRONOGICAL _ 

C  ORDER 

C  OUTPUT!  MEAN  ANO  STAND  DEV  3r  LAG  ANGLE  _ 

C  ALL  ANGLES  ARE  IN  UNITS  OF  RADIAN  *••*• 

C  TAU!  DELAY  IN  SECONDS  .  _ _ _ 

C  ALPHA!  PARAMETER  VECTOR 

C  ELERR!  MEAN  EL  LAG  ANGLE  (I.E.  TARGET  ANGLE-BARREL  ANGLE) _ 

C  AZERRI  MEAN  AZ  LAS  ANGLE 

C  ELSOI  STANOARO  DEVIATION  OF  ELEVATION  LAG  ANGLE _  _  _ 

C  AZSOl  STANOARO  DEVIATION  OF  AZ  LAG  ANGLE 

C  ELTRI  MEAN  EL  TRACER  ERROR  (TARGET  ANGLE-TRACER  ENDING  ANGLE) _ 

C  AZTRI  MEAN  AZ  TRACER  ERROR 

C  ELBARl  MEAN  EL  DARREL  ANGLE  _ 


C  OEL  >  TIME  STEP  USED  IN  THE  INTEGRATION  ROUTINE 

C  TS(  I)  !  STARTING  TIME  OF  I-TH  3 LANKING  DURATION  _  _  _ 

C  TE( I) I  ENDING  TIME  OF  I-TH  BLANKING  DURATION 

C  UO (1)1  INITIAL  GUESS  OF  EL  LAG  ANG.E  _  _ 

C  Y102I  INITIAL  GUESS  OF  AZ  LAG  ANGLE 

C  UELI  EL  CONTROL  _ 


C  UAZI  AZ  CONTROL 

C  C0(1>  I  EL  RATE  CONTROL  COEFF  _  _ _ 

C  C0(  2)  I  AZ  RATE  CONTROL  COEFF 

C  Kll  NO  OF  POINTS  IN  HE  ENTIRE  TRAJECTORY  _ 

C  Kl  NO  OF  POINTS  AFTER  THE  CIRST  TRACER  ROUNO  IS  FIREO 

C  EL 00!  EL  ANGULAR  ACCELERATION  OF  TARGET _ _ 

C  AZOOI  AZ  ANGULAR  ACCELERATION  OF  TARGET 

C  X3(l) !  EL  ANGULAR  VELOCITY  OF  TARGET  _  _ 

C  X3(  2)  I  AZ  ANGULAR  VELOCITY  OF  URGE! 

C  x«.l  ESTIMATION  ERPDR  OF  ANGULAR  VELOCITY  OF  TARGET _ _ 

C  ELI  EL  ANGULAR  “OSITION  OF  TARGET 

C  H(l) I  MODEL  PREOICTEO  LAG  ANGLE _ 

CH(2)I  MOOEL  PREDICTED  TRACER  ERROR 

C  Pit, 1M  VARIANCE  0 *  PPEDICTSD  LAG  ANGLE  _  _ 

C  P(2,2) l  VARIANCE  0T  PREOICTEO  TRACER  ERROR 
C  TO!  THE  INITIAL  FIRING  TIME 
C 

REAO*,K1»TO» IpT» IBL 


PRINT  3,K1»T0»I9T 

3  FORMAT  ( LM1, "NO  OF  PTS*  ",  I ,,  2X,"INIT  TIME*  *  ,512.9// IX, "REAO  EVERY 

C", 12."  POINT"/) 

IFU9L.GT.D)  READ*,  (TS(<).  TEIK),K»1,  IBL) 

PRINT  11, IBL 

11 _ FORMAT (lXj  I5j  IX,  “BLANKING  INTERVALS  ARE  “/> _ 

IFlIBL.CE.i)  PRINT  i,,ITS(0,TF(<0,K«l,lBL) 

%  FORMAT  <5I1X,"(",  F9.2,~,",F9,2,")“)I 
KT-T0/0.03 

K»K1-KT  _  _ 

T«TO 

NOIMM 
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}■ 


I  PRINT *20 /T^T 
C0(l>*1.34 
C0<2)=1.2« 

N01»N0IM-1 
N02*N0IM-2 

_ NAA=NOlH/2-l _ 

isr=i 

DEL=0 . 03* I°T 
ELSO*.Q05  **1.5 
AZSD*. 005**1. 5 
21*0. 

_ 21*JU _ 

ISET*0 
U£L=0. 

UAZ=  Q. 

PRINT  7 

7  FORMAT  (/1H  ,2X  ,  -  T I  ME**,  IX,  -  EL  VE.**,9K  ,**EL  ERR  ",  5  X,  ME  LSD”  *  6X, 

1**EL  CTR",6  X,”A?  VEL^ifeKt^MZERR^tfeXf^AZS CT  R** _ 

2,6X,~EL  TR'*,6X,  '*  A  7  TR'V) 

REAO  < 1 1  2)  <TfOlMi,AZtAZH,  AZQO,EL, ELO, ELOO.AZMN, <  , ZS  O  ,  S,  T*  1, KT 
C  > 

OO  5  1=1, K 

REAO(l,2) T,1UM1,AZ,AZD,AZ11,EL,EL3,ELDD, AZMN,X,ZSD,S 
_ IFCEOF(l) ) 1.1 _ 

1  IF<MOOa-l*IPT).N£ .  1)  SO  TO  5 
IH* ( 1-1 ) / I°T*1 
T*TQ ♦ < IN-1 1 *OEL 
TAO=T.5 

I F ( OUM1 .LE . 2677.  >  TAU*0UM1/ < 930. - . 19*OUM II 

_ HT  AU*T  AU/OEL _ _ _ 

IF ( 1ST • EQ. 1)  DTAU*TAU 
1ST* 1ST  *1  __  _  _ 

RA*NAA/TAU 

A2  =  1#-(TAU-9TAJI  /DEL  _ 

OT  AU*TAU 

_ X3( 11 *£LP _ 

X3(2)*AZO 

IFUIM-ll  .SE.MTAUI  03  0  9  _ 

TAUR=TAU 

C  IFCOUMl.LE. 44011  T AUR*T AU* ANAX1 ( 0 . 6, OUM1/503 0 . 1  _ 

Y10(1) *-TAUR*ELl *• 00 1* < 5.  2 * TAUR* . 465 *TAUR**2 1 *C OS< EL*0. 0 5 1 

C  Y1Q2»-0.025*SIjM (l.tAZO) _ 

Y102*-TAUR*A71 

IF(IM.NE.I)  GO  TO  13__  _ _ _ _ 

EL9=EL 

A  ZO  =AZ  _ _ 

E10*Y10<1> 

_ £?.&!  UJL2 _ 

GO  TO  10 

9  ISET-ISETH  _ 

IFaSET.NE.l)  S3  TO  10 

Zl*EL~ (EL  0-E10 )  ♦•90lM5.2*rAU*.486*TAU**2)*COS(  ELG  ♦0J.05I _ 

Z2=  AZ- ( AZ  0-E20I 

10  _ C  ALL  09SEL6 (ELERRt  ELSO) _ 

ELBAR*EL-EL£RR 

IFMIM-l)  .LE.MTAUI  Y10(2)«ri02*:OSI-LeAR) 

CALL  03SAZMA7ERR,  AZSO.EL34R) 

IF((MOO(IH-1,I»?IYT).E3.0I  .3R.IIH.E3.il)  PRINT  6,T  ,  X3I11,  ELERR 
1,ELSD,UEL«X3(2>, AZERR, AZSO,  UAZ , ELTR, AZTR 

A  FORMAT  <  11 G 12.  5) _ 

9  CONTINUE 

2  FORMAT (12G12.5)  _  _ 

STOP 

CHO 

SUBROUTINE  13SE.6(ELE»R*iLS0> 

COMHON/S/C3 (21  ,1EL ,<<,N3,  YiO (2>, X30( 21 , EL. EDO .A  ZOO ,M,RA, A 2 


o  o  o 


A,N01,N02,N4A,UfUA? ,  ELT*,  4 Z  TR,  ISET,  Z1 ,12,  TAU,  TS  ( 15  >  ,  T£(i$) ,T,IBL 
OINENSION  WC4)  fPCVt  4>,PL<4,6> ,P2(4,4>,  ALPHA(7)  ,  A(  4,  4) 

1,6(16)  tEBI16>,E6Wll6),EA(4,4>,EAlMTC4,4>,F<4)  ,CQC (4,4> , C< 4) 

2* X3I 1000) t X41 1 99  0) fEDH(1033l  , ED3H( 10  OG)  ,  THET (130  0)  , ALP5OU000) 
EQUIVALENCE  IA(1 ,1)  ,3(1) )  f  (EA(1,  1)  trB(l)  )•  (E  AINT(1  ,  1)  ,EMM ( 1) ) 

DATA  NT/1*/ _ 

DATA  ALPHA/1.5471,  *  017491,  • 024433, .42316, .22 446E-7, ♦17975E-3, 

A  .17302E-3/ _ _ 

INITIALIZATION  _ 

IFt(KK-ll.GE.H)  GO  TO  6 _ 

IF (KK.GT ft 1)  SO  TO  5 

N1*N0#*2  _ _ _  _ 

NOIH-NO 

ALPl*ALFlO*ALP-»m> _ _ _ _ _  __ 

ALP2-ALPHAI?) 

ALP3»ALPHA(3) _ 

ALP4*ALPH A( 4) 

ALP5*ALPHA(3)  _  _ 

ALP6*ALPHA (5) 

ALP7*ALPHA(7) _ _  _ _ 

ALP50(1>»AL®5* 

ACC*0. _ 

IFLA6»0 

SCAL*Cmi**2/a£L _ 1 _ _ _ 

00  1  1*1, NO 

00  i  J*  1,N3 _ _ 

P(I,J>*0. 

CQCJ  I, J) *6  » _ 

1  A  1 1, J)*0 » 

00  11  IslfNl  _  _ _  _  _ 

N(l)tO, 

0(I)*0.  _  _  .  _ _ 

F ( I) *0. 

ii  P(i,i)*o.oqio _ 

P(ltl)*0, 0030256279 

P  (2,21  *0. 0030  333677 _ _ _ _ 

X3(ll*X30(ll 

X4(  1)  *0.  _ _ _ _ _ 

EON (1) *X3 (1> 

_ EOON (1) *0  » _ 

S1*0# 

5  N  ( 1)  *Y13(  1)  _ _ _  _ 

6  IF(ISET  .EQ.l)  4(2>*Z1 

IF(  IBL.LT  •  1 )  53  T3  H  _  _  _ 

ALP5*ALPHA(5) 

IS»IFLAG»1 _ 

00  12  I*IS, I3L 

IF(T ,6E. TS (I) . AM 0. T •LT.fECI))  GD  TD  15 
*TT*ANlNt(1.5,A:C/3.) 

rf(T .GE.TE(I) .090*1 #LT •(T£(I)+ATT))  GO  TO  16 
XF(7.GE.(7Eim4TT>l  GO  TO  21 

_ ACC*0. _ 

IF  IT .LT «TS ( IS H  GO  TO  13 
GO  TO  12 

L6  ALP1*ALP10^(AL»*»AI1)-AL»U)  M 1  .•£<*( -C  ,4  3*  (T-TEC II  II) 

ALP*»-0. 0001*  (1* -EXP  (-0.43MT-TEC  III  II 
GO  TO  16 

15 _ ACC*  ACC  ♦PEL _ 

XFLAG*I-1 

ALPl*ALPMA(U*E<P(-0.0»5S*(T-TS(im 
ALP3*ALPHA(  1>*S<Pl-3.52MT  *TS  ( I)  I) 
ALP5*O.OOOlMl.-EXP(-9.1?MT-TS(im  > 

ALP10*ALP1 
*9  IP  46 . 
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AlPl«AlPHA<11 
ALP3*ALPHAn» 

ALP5*ALPHA(5> 

GO  TO  18 
12  CONTINUE 

18  CONTINUE _ _ 

ALP50<KK)*ALP5  f 

ARG»-0EL*AL*>1 

IFURG.GT.-200D  S1*EXP<AR;) 

TMET|«)=EL 
K1*KK*1 

_ K2«KK-1 _ 

C 

C  COMPUTE  TARGET  VELOSITV  ANO  ESTIMATION  ERROR 
C 

X3CKi>»X3<>«)*n0*0£l 
IFIALP1.EQ.3.)  SO  TO  4 

_ X4(K1)  *Sl*X4<KtQ  »EOQMl.-Sl)/Alal _ 

GO  TO  3 

4  X4(K1)  *X4  <«>  ♦  E*)0*DSL  _  _ 

3  CONTINUE 

6 OH ( KK) *X3(<K)-X4<K<)  _  _ 

IFCKK.GE.2)  EOOH(KKI*(EON(  <K)  •EOMI K? ))  /  DEL 

_ X70<1)»X31K1I _ 

C0R*C<MU*ALP2 
C R4*C0 ( II  #AUp4 
IF!  IKK-1)  »LE«MJ  SO  TO  150 

All  *1, +0,001* (5*  2*T4U*3«  43 G^T AU*  *2 ) • SIN ( THET  CKK^M)  ) 

A(1,1)**C0? 

A  <1#  2)  *«»C0(  1>  *ALP3 _ 

A<  2t  NOl)  *AC  It  1)*AL1*A2 
AC2tNO)*A(lt2)*All*A2 

00  2  X*lf N02 

Jl*I+2  _ _ __  _ 

A ( Jl« I) *RA 

A  <  Jl«  Jl) »*RA _ 

2  CONTINUE 

CALL  DSCRTfNOtSt  0EL«E8«EBIXTf 5)  _ _ _ 

CR3*CR4*AL1*A? 

SCAll*SCAl*UU*A2>**2  _ 

C 

C  COMPUTE  MEAN  TRACKING  ERROR  lUgj>  LAG  ANGLE) _ 

C 

U*AIP2»HU)  »AlQ1»rf(2l+ALP4»(X3<KK»>X4UKD _ 

00  110  I* 1»N0 

00  120  J*  1* NO  _  _ _ _ _ 

0(I>*0II)+EACIf J)*W|  J) 

120  _ CONTINUE _ 

110  CONTINUE 

F(1I»<1.-CR4) *<1<KKI>CR4*X4IKK) 

F  ( 2)  »X3(  KK)  *CR3*  U4  0CK-N)  -<3  <  KK-MI  l  +  0  •  00  1*  (1,-A2)  *  (  5.2*  0. 972'T AU 

ll*C0S<THETC«O<-MH  _  _ 

00  130  I* ItNO 

_ DOL 140  J«l»2 _ 

OlliiOd)  fEAINM  I,  J>*F<  J) 

140  CONTINUE 
NIDOtl) 

OII)«0*  _  _ 

130  CONTINUE 

C _ 

C  COMPUTE  COVARIANCE  MATRIX 

C  _ 

CQC(i»  ll*(Al»5+AL»6*ABS<E^4KK)>  ♦  4L* 7* A9S < E0DH(  KK>  >  > 

1  *SCAL 

COC  (2,2)*CALP«"M  K< -M>  +41*3* A3S  <EOM<<K-M)  )+AL  P7* ABS ( EOOMi KK-M) 

1>  PSCAL1 


CALL  NULMElIvr.CSC.NfVU,  >1,101 
CALL  MUlT(?1,P,>I0,N1,P2,13) 

DO  22u  1*1, MD 
00  223  J*1,M0 
P(I,  J)*P1II,J)»>2(I,  J» 

220  CONTINUE _ 

150  CONTINUE 

ELERR*M(1)  _ 

ELSD*SQRT  (>(1,1)  > 

ELTR*M(2) 

RETURN 

_ ENO _ 

SUBROUTINE  3BSA»6( AZERR, AZSD, ELJ  » 

COMHON/S/CO(2I  ,  3EL  »  «,  NO,  T 10  ( 2) ,  X3  0I2) , EL, EL  00, AZOO.M.RA, 

1  A2,N01,ND2,NA1,UEL,U,ELTR,AZTR,  ISET  ,  Zl,  Z2,  TAU,  TS(  151  ,TE  (15)  ,T,IBL 

DIMENSION  MIL)  ,>(4,4)  ,»1(4,4)  ,P2  14,1,)  ,ALPHA(7),  A(4, 4)  _ 

1,8< 16) ,£B (15) ,£T IVT ( 16) »E4  <  4 ,4 ) , EA IVT<4, 4) , F<4)  ,CQC (4,4) , 0(4) 

_  2,X3(1000)  ,X 4(1)3  0)  ,  EOH ( 1 001 )  ,  ED3H (  10  00)  ,  ALP50tlQ0C) _ 

EQUIVALENCE  ( »  ( 1 , 1 )  ,  3(  1) )  ,  (  EA  ( l,  1)  ,  EB(  1)  ),(EAINT(1,1),EBIM(1)) 

DATA  MT/1./  _ 

OAT  A  ALPHA/5. 5TR4,.11«94,.  17773, 1.0J53,. 2528 6E -5,.  19766E-3,. 757 «5E 

A- 3/  _  _ _ 

IFUKK-l)  .GE.M)  GO  TO  6 

IF(KK.GT.l)  SO  TO  S _ 

N1*N0**2 

C 

C  INITIALIZATION 

C  _ 

00  1  1*1, NO 

_  00  1  J*l, NO _ 

P(I,  J)*0. 

COC  (If  J)  =  0  •  _ 

I  A(I,JJ*0. 

00  11  1*1  »**> 

0(l)s0, 

F  < 1 1 *0 • 

II  P (It  I! *0*0190 

P<1, 1>*(0« 3)5918 3*C3S(£13>> 

P(2t2)*(0#1171587*CO$(El5) > **2 
X3<1>*X30(2) 

— _ 

EDH(  1)  »X3<  1) 

EOOHfllaO. 

Sl»9. 

ALP1*ALP10  * AL*HA  ( 1 ) 

AlP2*AtPHA«2) 

_ #tP3*AtPHAL(? _ 

AIP<»*AIPHAU) 

ALP5*ALPHA(5I 

AtP6*ALPHA(5l 

AlP7»ALPHAt7) 

ALP50( 1) *  AL°5 

_ _A  CC«0. _ 

iPLA6»0 

C 

C  CONPOTE  AND  STO*E  STATES  X3  A>ID  Xii 

C  _ _ 

5  M(l)«riO(2) 

6  _ IMISET.EQ.l)  4( 2)»Z2»:0S(£LG) _ 

CONTINUE 

Ca«COS(ELG) 

tP(IBL.LT.l)  GO  TO  Id 

ALP5*ALPHA(3> 

tS«IflAG*l 

09.12  l*XSt IBL 
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i  w  o  O  O  OOO  ooo 


IF(T.Gt*TStI>  •  0*  T  .UT.rECll)  SO  TO  15 

ATT*AMINl(l*«$f  A3C/3.) 

IF<T.Gt,Tnr».4ND#T.lT.<r-(IMArTn  GO  TO  16 
iF<r.G£«<r£<riMm)  go  ro  21 
ACC«0 • 

_ lF(T.LT»TStI$n  GO  TO  IS _ 

GO  TO  12 

16  ALPMALPIQMALMAU)  -AL^IOI  •  < 1 . - EXP< -0 . 4 3*  (T -TECI)  I  )  ) 
ALP5*-0.0QlMl.-EXP<-O*43Mr-TE<I>>)  > 

GO  TO  ia 

15  ACC=ACC*DEL 

_ |  FLAG»I-1 _ 

ALP1* ALPHA  Cl) *E<P<- t .075 5* ( T -T S( I) ) I 
ALP2=ALPHA(2)  *  EX  p  <  -  3 . 1 2M  T  •  T  S  (II  I) 

ALPJ*AL  P-<  AM)  *UP(-0,*>*tT-TSll>  )) 
alp5*0.ooimi.-exm-o*i2m  r*TS(nm 
ALP10=ALP1 

GO  TO  16 _ 

21  ALPlsALPHA(l) 

ALP2*ALPHA(2) 

ALP3  =  ALPHAM) 

ALP5-ALpH A( 5> 

GO  TO  15 

12 _ CONTINUE _ 

16  CONTINUE 
ALP50<<<MALP5 
ARG=-OEL*Al»lM3 

IF ( A9G  «GT  *-200 • )  Sl*EXo<ARD  _  _ 

Kl=«*l 

_ K2*KK>1 _ 

X3<K1MX3<«>  *A*DD*OEL 

IF (  ALPl »£0 )  jO  TO  4  _ _ 

X4<<1MS1MX4< «>  *AZ00*3ED 
GO  TO  3 

X4(K1MX4<K<>  ♦A^OP^OEL 

CONTINUE _ 

COMPUTE  ANO  STORE  ESTIMATED  TARGET  XEL03ITY  AND  ACCELERATION 
EDH<K1)=X3<<1) -XMKD 

IF<<2#Ge.  1)  Ep-H<«M(E3H<«)-E3H<<2)  )/OEL 

_ X30<2)  *X3  <  Kl) _ 

IF<|«-1)  .LE.M)  30  TO  153 

OO  2  1*1 « N02  _  __ 

J1 *1*2 
A  <  J  It  IMRA 
AU1#  J1M-RA 

CONTX_NUE  _ 

THE0O* (ELG-7EL3) /DEL 
CB*COS  <ELG> 

Tb* -TMEB0*T  AN  <  EL  G ) 

SCAL* <C 0<  2)  *C?I  *M/3EL 
C0R*C0 <  2)  *  ,lp?*:b 

_  A<lt  n»-C0?»TB _ 

A<lt2)*-:0<2)*ALP3*C3 
A<?f NOl)*-C)R*A? 

A(?tNO)*A<lt2)*A2 
A  <  2  *  2MTQ 

COMPUTE  TRA  NS  IT  TON  MATRIX  FA  AXP  ITS  INTEGRAL  EA INT _ 

CALL  OSCRT<NO,3,OELtE6tE3I><r,5) 

CR4  *Cfl<2)*ALP4#CB 

COMPUTE  MEAN  TRACKING  ERROR 
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U*ALP2*W(  1)  ♦  AL91*N<2)+ALP4MX3<<K>-<4<KK>) 

DO  110  I *  1  *  NO 
00  120  U*l*NO 
0<I)*0<I>+EAU,  J)*M< Jl 
120  CONTINUE 

11  0 _ CONTINUE _ 

FC1)*<CB-CR4)*XT<KK)  ♦CR!**X5dKK> 

F  t  2)  *X3(  KK)  •CBORUMXUfKX-N)  -X3<  KK-N) J  *A2 
00  130  1* l»NO 
00  140  J* It? 

OII)*DII) ♦EAINM I, J) ♦FIJI 

140 _ CONT INUE _ 

WU>*0<I> 

D  ( I) *  0  • 

130  CONTINUE 
C 

C  COMPUTE  COVARIANCE  MATRIX 

5 

COCTI  ,lT=TA  L^5  *AL  3  6#AB5CEDH(KKll  ♦  AL3  7* ABSTeFohTkK) ) ) 

1  *SCAL 

CQC (2*2)*(ALP5T(K<-N)*ALP6*ABS(;0H(<K-M) ) 

1  ♦ALP7*A&S(EPOH (KK-M) )  I  * S CAL* A2 *• 2 
CALL  MULT (SAINT,  CQC,U0,N1,  *1,101 

_ C A LL  MULTiEA , p, M P, Nlt P2, 13  > _ 

00"  220  1*1  *  NO 
00  220  J*  1 1  NO 
PCI*  J)*P1(I,JM*32(I,  J) 

220  CONTINUE 
150  CONTINUE 
OELG*ELG 


AZERR*H(l>/09 

AZSD*S0RT (Pll»l))/C9 

AZTR*W(2)/C3 

RETURN 

ENO 

SUBROUTINE  NULH  E  •  F,  L,  LI .  H.  NR ) 


01 ME NS I  ON  E(l) t c ( 1 ) #3(16)* •!( 1) 

00  10  1*1, L  _  _  __  _ 

11*1 

DO  10  K*1,L  _  _  __  __  _ 

TEMP*0. 

00  5  J*I ,  L 1 « L 


TEMP*T  EMP ♦E(J)*e(II) 

5  II*II*1  _  ___ 

<K*  (K-l  >  * L ♦> I 
HCKK)*TEMP 
10  G(KK)*TEMP 

IF(MR.EQ.l) return 


00  20  1*1 fL 

00  20  K*It l  _ 

T  EMP*0 • 

II*K  _  __ 

00  15  JaltLltt 

TEMP* TEMP+GC  J)  *S  (I  I) _ 

15  II*I I+L 

KK*(K-ll*l*t  _ 

20  M«KK)*TEMP 

L2*L-1 

00  30  1*1 «L2 

_ L3*XM _ 

00  30  J*L3tL 

Kl*(I-l)*l*J  _ 

K2*  ( J-ll *L*I 
30  H( Kll *HC  K2) 

ENO 

$U}ROyTINE  OSCRrtNOlM,AvoeL«EA*EAlNr«Nn 
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DIMENSION  A(i)r-A(l) ,£410 <1> ,CDEF< JO) 

C  SETS  EA*2X9(4*OEUtEAHr*INTSGRAL  EA  0  TO  DEL 

NOIHi  =  NOIMM 
NN=NOIM*NOIN 
NTM1*NT-1 

_ COEF<NT)«l> _ 

00  10  I=lt  NTM1 
I I *NT-I 

10  COEF(III*OEL*C3'F< II*l)/-l3AT< I) 

C  NT  MUST  3E  AT  LEAST  3 

CALL  DIAGmn,-AINT,A#:D£-  (1)  t:0EFC2)  I 

_  J>0  60  L>3tNT _ 

CALL  MULT(Af£AISTf NOIMt^N, - A,l) 

IF(L.£Q.NT)30  n  70 

60  CALL  OIAGCNOIM, EAINT, EAti • J#COEF(L> I 
70  00  80  ri=ltNN, ^3IM1 
EACH)  =EA(It)  *1*0 

_ 80  CONTINUE _ 

END 

SUBROUTINE  OT A3( NT IM f A , 3 f S 1» C2) 

DIMENSION  A ( 1) *3(11 
NDIH1*N0IM+1 
NN*ND IM^NDI 1 

_ N  MjsNDIM-1 _ 

21*1 

IF (Cl • £Q« 1 • 3 >  33  TO  10 
00  5  J*i»NN,MOM 
K* J+NMl 
00  4  I*J,K 

_ 4  A(I)*C1*9(H _ 

A ( II) * A( II) ♦C? 

5  1I*II*NQIM1 
RETURN 

10  DO  7  J=1,NN,N0I1 
K* J+NM1 

_  DO  6  I*  J»  K _ 

6  A( I) *B ( I) 

A  < II) *A  < II) +C2 

7  II=II*NQIM1 

RETURN  _ 

ENO 


i 
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